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Search for the Invisible Decay of Neutrons with KamLAND 
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The Kamioka Liquid scintillator Anti-Neutrino Detector (KamLAND) is used in a search for single neutron or 
two neutron intra-nuclear disappearance that would produce holes in the s-shell energy level of 12 C nuclei. Such 
holes could be created as a result of nucleon decay into invisible modes (inv), e.g. n —> 2>v or nn —> 2v. The 
de-excitation of the corresponding daughter nucleus results in a sequence of space and time correlated events 
observable in the liquid scintillator detector. We report on new limits for one- and two-neutron disappearance: 
r(n — > inv) > 5.8 x 10 29 years and r(nn —> inv) > 1.4 x 10 3 ° years at 90% CL. These results represent an 
improvement of factors of and > 10 4 over previous experiments. 

PACS numbers: 13.30.Ce, ll.30.Fs, 14.20.Dh, 29.40.Mc 


Baryon number violation is an important signal of physics 
beyond the Standard Model. Grand unified theories suggest 
that processes such as p —► e + 7r° and n —► K°D are important 
but suppressed by powers of the grand unification scale. How¬ 
ever, baryon number violation is not limited to grand unified 
theories. In fact, baryon number is an “accidental” symmetry 
arising from the pattern of particles in, and renormalizability 
of, the Standard Model. New particles or new physics such as 
supersymmetry or extra dimensions might also be the cause 
of baryon number violation. Moreover, the suppression may 
only be at the TeV scale, instead of the much higher grand- 
unified scale. 

The Particle Data Group lists more than 70 possi¬ 
ble modes of nucleon decay which conserve electric charge, 
energy-momentum, and angular momentum. Experimental 
lifetime limits of 1 O' 10 years for all but a few of these modes 
have been obtained. The decay modes with the poorest lim¬ 


its are so called “invisible” modes, such as n —> z/s and 
nn —> i/s. Invisible modes are dominant in some models (2l- 

Invisible decay modes involving bound neutrons are poten¬ 
tially detectable even though no energetic charged particles 
are produced directly in the decay. If we assume that the un¬ 
observed decay products carry away most of the neutron rest- 
mass energy, then the residual nucleus, k ~) X or A j?X, has a 
hole or holes in the previously occupied shell. Particles emit¬ 
ted in the nuclear de-excitation of the daughter are the experi¬ 
mental signature for the process. Detectability is independent 
of the specifics of the process as long as the rest-mass energy 
is carried away by the undetected particles. 

The SNO Collaboration established the best single neutron 
disappearance limit of r(n —> inv) > 1.9 x 10 29 years 
(90% CL) by searching for the de-excitation 7 rays follow¬ 
ing neutron disappearance in le O 01]. The best limit for 
two-neutron disappearance, r(nn —> inv) > 4.9 x 10 25 y 
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TABLE I: Branching ratios for the n C* ( 10 C*) 0 de-excitation 
modes after neutron (two neutron) disappearance from the s’i/2-state 
in 12 C and the experimental signature (number of time-correlated 
hits) for observation of these modes in KamLAND. 


Decay 

mode 

Daughter 

(decay, T 1/2 , and Q ec) 

Branching 

ratio 

% 

Exp. 

sign. 

(hits) 

(nl) xl C(n) 

10 C gs (/? + ; 19.3 s, 3.65 MeV) 

3.0 

3 

(n2) n C(n, 7 ) 

10 C gs (/? + ; 19.3 s, 3.65 MeV) 

2.8 

3 

(nnl) 10 C(n) 

9 C(/3+, 0.127 s, 16.5 MeV) 

6.2 

3 

(nn2) 10 C(n,p) 

8 B(/3+q, 0.77 s, 18 MeV) 

6.0 

3 


(90% CL), was set by the Borexino Collaboration 01 by 
searching for possible decays of unstable nuclides resulting 
from nn disappearance in 12 C, 13 C, and le O. 

KamLAND, currently the world’s largest low-background 
liquid scintillator (LS) detector (the LS composition is 
CH 1 . 97 ), has a low energy threshold (< 1 MeV) and good 
energy and spatial resolution, making it well suited for neu¬ 
tron disappearance searches. The LS is contained in a 13-m- 
diameter transparent nylon-based balloon suspended in non¬ 
scintillating purified mineral oil. The scintillator is viewed by 
1325 17-inch Hamamatsu and 554 slower 20-inch Hamamatsu 
photomultiplier tubes (PMTs) mounted on the inner surface 
of an 18-m-diameter spherical stainless-steel tank. The in¬ 
active mineral oil buffer serves as passive shielding against 
external backgrounds such as 208 T1 7 rays coming from the 
PMT’s glass and nearby rocks. A 3200 metric ton water- 
Cherenkov detector surrounds the stainless-steel sphere, act¬ 
ing as a cosmic-ray muon veto. The muon rate in the central 
detector is 0.34 Hz. Details of KamLAND can be found else¬ 
where a. 

In this study, we consider the disappearance of neutrons 
only from the fully occupied s-shell in 12 C leaving the daugh¬ 
ter nucleus in a highly excited state. This results in sub¬ 
sequent de-excitation with the emission of secondary parti¬ 
cles (p. n. d. a , 7) and the possible further decay of the resid¬ 
ual radioactive nuclei. The probabilities for the various de¬ 
excitation modes following one- and two-neutron disappear¬ 
ance in 12 C were estimated with a statistical nuclear model by 
Kamyshkov and Kolbe 0. Some of these modes, producing a 
sequence of time-and-space-correlated events, are detectable 
in KamLAND. Here we report on a search for the four de¬ 
cay modes listed in Table H] each producing a correlated triple 
signal in the detector. 

According to Ref. 0, neutrons produced in these modes 
have energies up to 45 MeV. In the LS, the first scintillation 
signal arises from neutron-proton elastic scattering, 4.4 MeV 
7 rays from the inelastic neutron scattering from 12 C, and de¬ 
tectable secondaries of other neutron interactions in the scin¬ 
tillator. Neutrons thermalize and capture, mostly on protons 
(^capture ~ 210 /ts), producing a second 2.2 MeV 7 ray sig- 


TABLE II: Top section: Selection criteria used in the neutron (two 
neutron) disappearance search. Bottom section: detection efficiency 
e for the decay modes nl, n2, nnl and nn2 given in TableQ] 


Quantity 

n disappearance 

nn disappearance 

Rl,2,3 

[m] 

5.0 

5.5 

Rxy3 

[m] 

>1.0 

>1.0 

AR12 

[m] 

2.0 

2.0 

AR13 

[m] 

0.8 

1.0 

AT12 

[ps] 

0.5 -1000 

0.5 -1000 

at 13 

[s] 

0.003-70 

0.003 - 6 

Ei 

[MeV] 

0.9 - 25 

0.9-40 

e 2 

[MeV] 

1.8 -2.6 

1.8 - 2.6 

e 3 

[MeV] 

1.5-3.8 

3.1-18.0 

Lil(nnl) 


0.430 ±0.027 

0.680 ±0.032 

e n2(nn2) 


0.651 ±0.033 

0.678 ±0.032 


nal. A third correlated signal comes from the /3 + decay of the 
residual nuclei ( 10 C, 9 C or 8 B). The correlated triple signal is 
a powerful tool for identifying these decays in the presence of 
backgrounds. 

The limits presented here are based on data collected be¬ 
tween March 5, 2002 and October 31, 2004. Initially, 
KamLAND had 17-inch PMTs only operational with 22% 
photo-cathode coverage and 7.3%/ \J L’(MeV) energy resolu¬ 
tion. After February 27, 2003 the photo-cathode coverage in¬ 
creased to 34% by including 20-inch PMTs, which improved 
the energy resolution to 6.2%/i/£l(MeV). The location of 
the event in the detector is reconstructed from the arrival times 
measured with the PMTs. Typical spatial resolution for the re¬ 
construction procedure at the detector center, for a point-like 
energy deposition, is ~ 20.6 cm (~ 24.0cm)/y / E(MeV) with 
20-inch (without 20-inch) PMT data. The energy of the event 
is reconstructed using the number of detected photo-electrons 
after correcting for position and PMT gain variations. The 
response of the detector is studied and monitored with peri¬ 
odic deployments of 7 ray ( 203 Hg, 68 Ge, 65 Zn, and 60 Co) and 
neutron sources (Am-Be) along the central vertical axis of the 
detector. Natural radio-isotopes remaining in the detector ma¬ 
terials such as 40 K and bismuth-polonium from U/Th chains 
are also used for calibrations. Because of the external 7 ray 
background, a radial fiducial volume cut of 5 m (5.5 m) for 
the n (nn) disappearance search is applied. This corresponds 
to a fiducial mass of 408.5 (543.7) metric tons, correspond¬ 
ing to 3.48 x 10 31 (4.63 x 10 31 ) s-shell neutrons in 12 C nuclei. 
The total livetime, after correction for muon-associated cuts 
to avoid backgrounds due to muon-spallation events, is 749.8 
days (751.4days) for the 5.0m (5.5 m) fiducial volume. 

The correlated triple event-selection cuts are summarized 
in Table [d] The indices refer to the first, second and third 
energy deposits (hits) in each candidate event. R is the radial 
position of the hit, Ryv is the distance to the vertical central 
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FIG. 1: n disappearance search: (a) correlated triple background 
events in the 150 s-1000 s off-time window for the 3 rd hit inside the 
5.5 m detector radius; (b-d) events inside the 5 m fiducial volume and 
outside of the 1 m cylinder cut around the central axis are compared 
to the expected signal for r = 4xl0 27 y, generated with the Monte- 
Carlo. Dashed lines indicate cuts used in the data analysis. 


FIG. 2: nn disappearance search: (a) correlated triple background 
events in the 10 s-1000 s off-time window for the 3 rd hit inside 6 m 
detector radius; (b-d) events inside the 5.5 m fiducial volume and 
outside of the 1 m cylinder cut around the central axis are compared 
to the expected signal for r = 10 28 y, generated with the Monte Carlo. 
Dashed lines indicate cuts used in the data analysis. 


axis, E is the reconstructed energy of the hit, AR and AT 
are the position and time differences between the hits. After 
applying all cuts to the dataset, the number of observed n (nn) 
disappearance candidate events is 1 ( 0 ), compatible with the 
expected background, as discussed next. 

No other process is expected to produce triple events with 
sufficient rates to result in backgrounds for the present search. 
Backgrounds to the n and nn disappearance search originate 
from doubly correlated signatures, such as those derived from 
reactor T7 e events or other processes such as neutrons pro¬ 
duced in 13 C(a,n) 16 0 followed by an uncorrelated coinci¬ 
dence within the selection cuts. The third uncorrelated hit 
is dominated by the tail of the 210 Bi /3-spectrum and 208 T1 
7 rays (see Fig. m and Fig. ES). The data analysis cuts were 
optimized to suppress the acidental coincidence background 
and to improve the signal/background ratio. For the case of 
the n disappearance search. Fig. ns shows the location of the 
accidental background events within the 5.5 m fiducial volume 
using a relaxed AR 13 < 2 m cut. Events inside the 5 m fiducial 
volume and outside of the 1 m cylinder cut around the central 
axis (to remove background from thermometers positioned at 
0m and ±5.5m) are plotted in Fig. 0b, c and d) to show the 
position of the selected cuts. Similarly, Fig. 0 illustrates the 


selection criteria for the nn disappearance search. To remove 
the background from muon induced neutrons and after-muon 
detector activity, a 2 ms veto of the entire detector volume 
following a muon traversal is applied for the case of the 3 rd 
hit. The /3-neutron emitter 9 Fi is produced by muon spallation 
and can give prompt and delayed signals. This background is 
suppressed by a 2 s veto within a 3 m radius cylinder around 
a well-reconstructed muon track. The entire volume is ve¬ 
toed for 2 s following energy depositions with more than 10 6 
photo-electrons above what is expected from a minimum ion¬ 
izing particle or if the muon is not well reconstructed. For n 
disappearance the accidental background is estimated with a 
150 s to 1000 s delayed-coincidence off-time window for the 
3 rd hit. In the nn disappearance case a 10 s to 1000 s win¬ 
dow is used. The number of expected accidental background 
events in the entire data set for the n disappearance search is 
0.82±0.26. For the nn disappearance search, the significantly 
greater energy of the 3 rd hit (see Fig. 01) allows selection of 
the energy threshold cut to be larger than the energy of the 
main source of background; the number of expected acciden¬ 
tals is 0.018±0.010 events. 

The efficiency of the position and energy cuts is calculated 
using a Monte-Carlo procedure based on the GEANT code 0] 
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FIG. 3: KamLAND Am-Be data compared to the GEANT detector 
simulation: (a) prompt energy signal; (b) diffusion of the 2.22 MeV 
delayed 7 ray from neutron capture on proton 


in Sec. II-III of Ref. IT2I0 are used, xu m = 3.07 (xi vm = 2 . 28 ) 
events at 90% CL, for neutron (two neutron) disappearance 
respectively. The corresponding lifetime limits at 90% CL are 
r(n —> inv) > 7.2 x 10 29 y and r(nn —> inv ) > 1.7 x 
10 3 °y. If the Feldman-Cousins procedure III is used instead 
of the simple Poisson limit, xu m = 3.53 events at 90% CL 
(xum = 2.41 events at 90% CL) and the corresponding lifetime 
limits are slightly reduced: r (n —► inv) > 6.3 x 10 29 y 
and r (nn —► inv) >1.6 x 10 3 ° y. 

The effect of the uncertainty in the theoretical branching 
ratios and uncertainty for the number of expected background 
events is estimated using the program POLE 0. Calcula¬ 
tions using POLE are consistent with the results derived us¬ 
ing the Feldman-Cousins procedure in the case of no sys¬ 
tematic uncertainties and give xu m = 3.82 events at 90% CL 
(xum= 2.75 events at 90% CL) when the systematic uncer¬ 
tainties are included; the lifetime limits are: 


set up to simulate the KamLAND detector response to the de¬ 
cay of the n C* ( 10 C*) nuclear system. The code includes 
new libraries 0 which give a more accurate description (com¬ 
pared to standard GEANT) of the (n, 7) and (n, n' 7) reac¬ 
tions. The decay schemes of 10 C and 8 B are well known, 
the 9 C decay scheme is less well studied and the most re¬ 
cent data B2D are utilized for our simulations. The observed 
non-linear dependence of the energy estimation upon particle 
energy and type (for 7 rays and the neutron Am-Be source) 
is reproduced with a combination of a correction for scin¬ 
tillator light output quenching 0 and a theoretical correc¬ 
tion for Cherenkov light emission for electrons and positrons 
(Fig- lit)- The energy and position are taken from the Monte 
Carlo simulation and then smeared according to the resolu¬ 
tion observed in calibration 7 ray and other data (Fig. 03). 
The reliability of neutron simulations in the LS was tested by 
comparison to simulations using the SCINFUL code 0 and 
a good agreement was found. The calculated efficiencies for 
the de-excitation modes in Table [flare given in Table ETI 
The systematic uncertainty of the predicted theoretical nu¬ 
clear de-excitation branching ratios as estimated in 0 is 
rather large, ±30%, while the experimental systematic uncer¬ 
tainties, such as the fiducial volume uncertainty and the un¬ 
certainty in neutron energy scale estimated by comparison of 
the Am-Be data and GEANT simulations (see also Table un 
and Ref. ( 5 )]), are significantly smaller. The total systematic 
uncertainty is ~30%. The lifetime limits r(n(nn) —> inv) 
for n and nn disappearance are calculated using 


r(n(nn) -> inv) > — - ^ e n(nn ) J B n(nn)i , (1) 


%lim 


i= 1,2 


where Nq is the number of s-shell neutrons (or n pairs), T 
is the detector livetime and e n ( nn )i is the detection efficiency 
(given in Tableim. /i n(nn)? is the branching ratio from Table[j] 
and xum is the largest number of events compatible with ob¬ 
servation at a given Confidence Level. If classical confidence 
intervals for Poisson processes with background (as described 


r(n -*• inv) > 5.8 x 10 29 y at 90% CL (2) 

r(nn —► inv) > 1.4 x 10 3O y at 90% CL (3) 

To summarize, about 838 metric ton-years (1119 metric 
ton-years) of KamLAND data have been analyzed in order to 
search for disappearance of a single neutron (a neutron pair) 
from the s shell in 12 C which create subsequent nuclear de¬ 
excitations leading to three time and space correlated events. 
The observed number of events is consistent with that ex¬ 
pected from the accidental coincidence background. This re¬ 
sults in new improved limits on single neutron and two neu¬ 
tron invisible decay lifetimes. This KamLAND result is a fac¬ 
tor of 3 better than the present best limit from the SNO col¬ 
laboration Qj for n disappearance, and more than 4 orders of 
magnitude better than the nn disappearance limit set by the 
Borexino collaboration 0] . 
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